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Abstract 

We study the stability of the system of the Euler equations in the 
neighborhood of the stationary solution associated with the quasi isobaric 
profile in a gravity field. This situation corresponds to a Rayleigh- Taylor 
type problem with a smooth base density profile which goes from to 
p a (of Atwood number A = 1) given by the ablation front model with a 
thermal conductivity exponent v > 1. This linear analysis leads to the 
study of the Rayleigh equation for the perturbation of the velocity at the 
frequency k: 

~^ P °^^ + k '^ po ^ ~ ^Po( x )]u = 0. 

We denote by the terms 'eigenmode and eigenvalue' a L 2 solution of the 
Rayleigh equation associated with a value of 7. Let Lq > be given. The 
quasi isobaric profile is po(x) = p a £(-£-), where £ = — £). We prove 

that there exists L m (k), such that, for all < Lq < L m , there exists an 
eigenmode u such that the unique associated eigenvalue 7 is in [01,02], 
ot\ > 0. Its limit when Lq goes to zero is \fgk. We obtain an expansion 
of 7 in terms of Lq as follows: 



7 = 



yfl + 2(r(l + I))-i(25£n)£ + O((fcio) min(1 '* ) ) 

We identify in this paper the expression of the next term of the expansion 
of 7 in powers of Lq . 

Using the existence of a maximum growth rate A and the existence of at 
least one eigenvalue belonging to ]4, A[ (thanks to a semiclassical analy- 
sis), we perform the nonlinear analysis of the incompressible Euler system 
of equations using the method introduced by Grenier. This generalizes the 
result of Guo and Hwang (which was obtained in the case po(x) > pi > 0) 
to the case where po — » when x — *• — 00 and ko(x) = P p °^ satisfy ko 

regular enough, bounded, and kop 2 bounded, which is the case in the 
model associated with the quasi-isobaric profile, according to v > i. 
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Statement of the problem and main result 



In this paper, we study a theoretical system of equations deduced of the fluid 
dynamics analysis of an ablation front model. Such models have been studied 
from a physical point of view by many authors (see H.J. Kull and S.I. Anisimov 
|12) , V. Goncharov, [B], P. Clavin and L. Masse JT7] ) . They can be considered as 
a generalization in the ablation case of the Rayleigh- Taylor instability, studied 
in the pioneering works of J.W. Strutt (Lord Rayleigh) TJ5] and G. Taylor [2D] . 
The Rayleigh equation models the Rayleigh- Taylor instability. It is obtained by 
considering the linearization of the incompressible 2d Euler equations around the 
solution (p (x), 0, 0,pa(x)) (density, velocity, pressure) with -p + po(x)g = 0. 
The system of equations write 

d t p + d x {pu) + d z { P v) = {) 

dt(pU) + d x { P U 2 + P) + d z ( P UV) = -pg 

dt(pv) + d x (puv) + d z { P v 2 + p) = o UJ 

d r U + d z V = 



Write p = p + a, U = t>i, V = v 2 , P = Po + p, the linearized system is 

(2) 



dx 

p (x)d t vi + d x p = -ag 
p (x)d t v 2 +d z p = 
d x vi + d z v 2 = 0. 



from which one deduces, using v\ = ue lkz , the partial differential equation 

d , , , d 



Introduce 



d^ PQ ^^ d t 2 ^ + k2 Po{x)d 2 2 u = gk 2 p' Q {x)u. 



T ( X , z , t) = Q(x, z, t) P{X ' Z ' t] - MX) 



p(x,z,t) ' p (x) 

the system (JlJ is equivalent to 

diT + U.VT = k Q (x)uT 

d t U + {u.V)U + TVQ + TQk (x)e 1 = (1 - T)g (3) 
divLT = 0. 

It is a consequence of the equality Tp^ 1 Vp = TVQ + k TQei + Tg and of 
d t T + U.VT = -p- 2 (d tP + UVp). 

The associated linearized system in the neighborhood oi U = 0, T = 1, Q = 
is 

d t T = k (x)u 
divu = 

d t ~u + P Q 1 V{puQ)+fg = 0. 
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Assume that the perturbation is written as (the real part of) a normal mode 

e^e zkz u(x,kL ), 

where k is the wavelength of the transversal perturbation and 7 is the growth 
rate in time of this perturbation. We obtain the Rayleigh equation ([JJ (see C. 
Cherfils, P. A. Raviart and O.L. [3]): 

" Tx {po{x) lt + {k2po(x) ~ !^o(z)k(z) = 0. (4) 

We consider a family of density profiles po(x) such that po(x) — p where 
I/o is a characteristic length of the base solution. In one of the physical applica- 
tions, namely the case of the ICF, its magnitude is meters, hence allowing 
us to consider the limit Lq — > 0. 

We develop here a constructive method for the study of the modes associated 
with the Kull-Anisimov density profile (see B. Helffer and O.L. [H]). The Kull- 
Anisimov profile po is given by 

Po(x) = Pa^(-f-), (5) 

where the function £ is a non constant solution of 

e = r +1 (i-o, (6) 

v is called the thermal conduction index. 

Note that this equation on the density is NOT obtained from the incompress- 
ible Euler equations but from a compressible model with thermal conduction 
introduced by Kull and Anisimov [T2] and used for example in [TT] or in [T4] , 
The Kull-Anisimov profile satisfies lmi x ^ +oa po(x) = p a , where p a denotes the 
density of the ablated fluid, and the convergence is exponential, whereas 

lim^-oopo^) = 

and the convergence is rational ((— x)» po(x) — ► Co > when x — ► —00). The 
associated Atwood number is thus 1. Remark also that all non constant solu- 
tions of © differ from a translation. 

This case may be related to the case of the water waves (the density of air 
being much smaller than the density of water) . It is thus a limit case in all the 
theoretical set-up used for the study of Euler equations for fluids of different 
densities. 

Note that, in this case, the self adjoint operator associated with the equation 
(fj| is not coercive in H 1 (R). The methods of [3J, [IH] cannot be used directly. 
Moreover, the properties of po do not allow us to apply [3], because it relies on 
Po{x) > pi > 0. 

However, consider ko(x) — introduced in the abstract. In our case, it is 

equal to Lq £ v (1 — £), hence it is a continuous bounded function which admits 

a maximum L~h, and, for v > |, kop 2 is bounded. These properties are (for 
a more general profile) what is needed to obtain the nonlinear result. 
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Remarks Define the function r(t, e) through: 

-M-'-Wil a-b) = ~ f + eh-^rfjt, e). (7) 

£ £ £ VI 

There exists to > and £q > such that r(t, s) is bounded for t > to, < e < £q, 
and has a C°° expansion in e,e~ . Define S through 



e~S'(t, e) = e^^tHt ~ 7Z> lim t ^+oo5(t, e) = 0. 



-i^(-j) 1 

£(-!) ^ 

We have the identity 



£(-J)(^exp( £ ^ )£ )) = l (8) 

which implies that there exists a function r bounded for t >to and e < £o such 
that 

exp(— vs~S(t, e)) = 1 + £~t~~r{t, e). 
Let tt(y) = u(Lqd). The Rayleigh equation rewrites 

- + ( £ ^(y) ~ ^'(y)My) = o, (9) 

where e = fcLo and A = We will consider this equation from now on. 
We shall introduce two equivalent versions of this equation, which are: 

1. the system on (U + , V + ) such that U + (y,e) — u(y,e)e ey and V + (given by 
the first equation of the system below), v(y,e) = V + (y£)e~ £y : 



dU+ 

dy ~ °^ /v ^+ t 5(y) 



^- £ (A + l)y + + £ (l-A 2 K(y)?7 + , 



e(l - X)U+ + jfaV. 

I dy 

2. if we introduce w = the system on (u, w) is 



= (A 2 - l)u - \w + (i + e^5'(t, e))«;. 



(10) 



(11) 



The first part of the main result of this paper was presented in [13] . and the 
case where £(y) — £(l)(y + 1) _ ~ for y > was solved in [3J. The case of the 
global ablation system was treated in p3] and is published [T5] , 
We finally recall that, if there exists a solution in L 2 (R) of ([9]), then A satisfies 
the inequality (see [T0]1^1 

A>max(l, £ ^-^ ). (12) 

The main result of the first part of this paper is 



1 It is a consequence of max(|) = r^py 
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Theorem 1 1. There exists e$ > 0, and C° > such that, for all e €]0,£o[ 
there exists X(e) € [^, |] such that the Rayleigh equation {P]) admits a 
bounded solution u for X = X(e), which corresponds to the eigenmode u 

and the eigenvalue j(k,e) = an d X(s) satisfies 



|A(e)-l| <C°e-. 



2. We have the estimate 

A(e) =1 
= 1 

with a = min(l, ~). 

3. We have the expansion 



2(2s)*(T(l 

2(f)^(r(i 



1 - A( £ ) = 2(i)^(B (0))- 1 [l + + 



v 

00 — -2s 



where Bq(0) = — 2 s^e 
below in Proposition^ 



ds = -2"*r(l 



This result is a result, for k fixed, in the limit Lq 

k 



and Co (1,0) is calculated 
0. It writes also, for k fixed 



and for L < ^ as 



l + 2(f)i(r(l + i))-i+o( £ i) 



(13) 



Note that, in this case, the order of magnitude of 7— \fgk~ is not in kLo as in [3], 
but the result of [TDJ, based on po — p a ^x>o G for all 8' > is pertinent. 

We have also a result for k going to infinity, which can be stated as 

Proposition 1 a) Any value X(e) such that {P]) has a L 2 non zero solution 



satisfies 



kg 



(my 



< A 2 , where A 2 = j^—. 



b) Any sequence k — > satisfies the following 



,, Me) T 



ay) 



It is proven in [ID] . 

Remark that formula (fl3]) and Proposition [I] are not in contradiction. They 
lead to two different stabilizing mechanisms induced by the transition region: 
one is a low frequency stabilization when Lq — > and the other one is a high 
frequency stabilizing mechanism when k — > +00. It is important to notice that 
Propositions[2]and[3]below allow us to construct an (exact) solution u(y, A(e), e) 
of the Rayleigh equation hence giving an unstable mode 



u{x, z, t) 



*'u( — ,kL ,X(e))eVm 
^0 
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solution of the linearized Euler equations. Moreover, from Proposition (TTJ), one 
has the following: There exists k > 1, A(e) and u(y) such that e — ULq, u 

solution of ©, 7(fc,e) = yfefe, | < 70,e) < A, ||w(?/)|| L 2 = 1, u(0) > 0. 
From the construction of this particular solution, we deduce a nonlinear result. 
For simplicity, in what follows, we will denote by y(k) the eigenvalue j(k,Lo). 
From u, one deduces a solution 

U = «i, Q 1 ,T 1 )e ikz +^ t ] = R[(u(aO, -l u '(x), - ^^(a,), ^ u(a; )) e ^+7«'] 



of the linearized system. We thus consider a function V N = (0, 0, l)(Loic) + 
Ef=i <^(x, y, t) satisfying (^(F) - S N+1 R N +\ V N (x, z, 0)-(0, 0, a, l)(i x) 
5U(x, z, 0). We also construct the solution V(x, y,t) of the Euler system such 
that Emod{V) = and V(x,z,0) = (0,0, 2a, l)(L ar) +SU(x,y,Q). Introduce 
finally V (x,y,t) = V(x,y,t) — V N (x, y, t). This procedure constructs a solu- 
tion of the nonlinear system. 
We have the 

Theorem 2 1. There exists two constants A and Co, depending only on 
the properties of the Euler system, on the stationary solution and on the 
solution u(x), such that, for all 9 < 1, for all t e]0, In S( j oA [, one has 

the control of the approximate solution V N in H 3 , namely 

,n -hi , n.-tN „N 1 1 ^ SACoe-rW 



|T iV -l|| Hs + ||^ v ||^ + ||Q iV -go||^ <C 



1 - SACoeMV* 



and the leading order term of the approximate solution is the solution of 
the linear system as follows 

P 7(*)* 

iir* - > mm^e** - Acic i _ SACoel(k)t 

2. There exists No such that for any N > iVo, the function V d is well defined 
for t < ^p^y In | and satisfies the inequality 

7(k) o 

3. We have the inequality, for eo < Co A 



This paper is organized as follow. The sections 1, 2, 3 study the linear system 
and identify the behavior of the growth rate 7(fe) when Lq — > by constructing 



G 



the Evans function, and Section 4 constructs an approximate solution of the 
nonlinear system of Euler equations. 

We identify in a first section the family of solutions of ^ which are bounded 
when y — > +00 and we extend such solutions, for (e, A) in a compact B, on 
{{;~ 1 (eR), +00 [, where R is a constant depending only on B (Proposition . In 
the second section, for all to > 0, we calculate a solution of © which is bounded 
on ] — 00, — zg] (Proposition [3]) . 

A solution it of © which is in L 2 (M.) goes to zero when y — ► +00 as well as 
when y — ► —00. Moreover, as po(x) is a C°° function on M, any solution u of 
© is also in C°°. 

Notice that lim £ ^o(— £ £,~ 1 (( £ R)~)) = 77^: from which one deduces that there 

exists to such that < to < |lim £ _>o( — e£ _1 ((£-R)^))- 

The regions ] — 00, — ^] and [£~ 1 ((£.R)i), +00 [ overlap and 



AevR levW 



Hence the solution u belongs to the family of solutions described in proposition 
[2] (of the form C*u+(y,£)) and belongs to the family of solutions described in 
proposition [3] (of the form C**U(—ey,e)), that is 

u(v) = C.u+(v,e),y>€- 1 ({ER)$) 
u{y) = C^U{-ey,e),y<- t -f 

From the continuity of u and of u', one deduces that, for all y± G [- 



(corresponding to t x = -ey± G 4^]), we have C >t u + (y_ L , e) = C M J7(tx,e), 



1 3 

C^u+(y ± ,e) = -C»eU'(t ± ,e) 
Introduce the Wronskian (where £ _1 has been added for normalization purposes) 



W(y) = e-\u + {y, s)jL(U(-ey, e)) - ^(u + (y, e))U(-ey, e)). 
ay ay 

It is zero at y± — —et±. Conversely, if A and e are chosen such that the 
Wronskian is zero (in particular at a point y± = — ^), the function 

, J C**U(-sy,e),y<y±_ 

U{y) = < r U(-ey ± ,e) r , > (14) 

is, thanks to the Cauchy-Lipschitz theorem, a solution of ([4]). Moreover, it 
belongs to L 2 (R) owing to the properties of u + and of U. 

In Section [31 we compute the function W. As U and u + are solutions of the 
Rayleigh equation, which rewrites 



— {u+(y e)) = -^1^!± + ( £ 2 _ gA ^M ' 
dy 2 + ' £(y) dy £(y) ' 



f (lOMw) = C(2/o)W(y ) for all y, y I I - r > I 



the function W is solution of = — T^rVV, which implies the equality 
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This Wronskian can be computed for y± £ [~i^7r: ~Wi7r) usm S the expressions 
obtained for U and u + . We prove that it admits a unique root for < e < £o 
and A in a fixed compact, and we identify the expansion of this root in e, hence 
proving Theorem [T] Precise estimates of this solution are given in Section 3. 
In Section 4, after proving a H s result on a general solution of the linear system 
(taking into account a mixing of modes), we calculate all the terms Vj of the 
expansion of the approximate solution, the perturbation of order 5 being an 
eigenmode with a growth rate 7 g]^, A[, where A 2 = maxfco (%) -g^ ■ 

1 Construction of the family of bounded solu- 
tions in the dense region. 

The system jTOj writes -^U + = eM (£_(y), X)U + . When y -> +00, the matrix 
converges exponentially towards Mo(l,A), which eigenvalues are and 2, of 
associated eigenvectors (1, A — 1) and (1, A + 1). 
It is classical that 

Lemma 1 There exists a unique solution (£/+, V+) of H10\) which limit at y — > 
+00 is (1, A — 1). Moreover, there exists £0 > such that this solutioi^ admits 
an analytic expansion in e for £(y) 6 [£q, 1[- 

The proof of this result is for example a consequence of Levinson [16] . 

The aim of this section is to express precisely the coefficients of this expansion 

when £(y) — > and to deduce that one can extend the expression obtained for 

We consider, in what follows, the change of variable 



( = IW (16) 



We prove in this section the 



Proposition 2 Let K be a compact set and A G K . There exists e$ > and 
R > such that, for < e < Sq, the family of solutions of 110\) which is bounded 
when y — > +00 is characterize^, for y such that £(y) > (sR)i , by 

f U + (y,e) = l+(^§^(A((,s) 

\ V + (y,e) = \-l + (l-X)(l-0(B((,s). 

The associated solution of HP is u + {y,e) — U+(y, e)e~ £y . 



2 It can also be shown that there exists a unique solution (U , V) such that (U, V)e 2ey 
(1,A + 1) 

3 a general solution is K +{U, V) where K+ is a constant 
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Proof of Proposition [2] We write the analytic expansion in e: 
U =l + J2 s j Uj,V = A - 1 + £1v j- 



We deduce, in particular, 



^ = M(i-e(y)) 

dy 



^ = (A 2 -i)(i-e(y)) 



hence assuming Ui,Vi — ► when £ — > 1 (which is equivalent to dividing the 
solution by its limit when £ — > 1) we get 



_i-Ai-r +1 
«i - — "err 

l-A 2 1-r 



I V1 

The following recurrence system for j > 1 holds: 



= !(„,- (A -1)**) 

%i = (A+i)(« i -(A-i)e« J ) 



Usual methods for asymptotic expansions lead to the estimates (which are not 
sufficient for the proof of Proposition [2]) 

MA j 

\uj(y)\ + \vj(y)\ < • 
so 

However, using the relation 1 — £ = |£pr, we obtain the following estimates: 
Lemma 2 Let £o > given. For all j > 1, introduce aj and bj, such that 

TTie functions aj and bj are bounded, analytic functions of £, /or £ £ [£o>l]- 
TTiey satisfy 

\aj(0\<AR j ,\bj(0\<AR j , (18) 
where R depends only on A. 

We prove Lemma [2] by recurrence. Assume that this relation is true for j. We 
have the relations 

f ^ = (l-A)(6,-(A-l)a J ) i ^| w 

I ^ = (1- A)(A + l)(6 j - (A - l)a,) | ^4 TT 

from which we deduce, using the limit at £ — * 1 



u i+ i (y) = (1 - A) / ^(j+D+2 
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and 



We thus deduce that S, u ^ +1 ^Vj + i (y) and £ v V +1 > +1 Uj+i (y) are bounded functions 
when £ g]0, 1]. Moreover, if we assume \bj\ < AR? and |a,-| < Affl , then 



drj 



K-+i|< ARi|l-A|(|A-l| + l)£ 

< AR?\1 - A||A + 1|(|A -11 + 1) £ — 4l 



We end up with 

\u j+1 \<\x-i\AR ^-^r:S^\ 

|, J+1 |<|A-l|^|A + l|i™^ 



IK 



As ^ < 1- G [0, 1], we get < |A- < 

Ai#|A - 1[ ]A + I] (|A ~ 1 ^0~ C(y)) ■ Consider 

i? A = (|A-l| + l)max(l,|A + l|). (19) 
The previous inequalities become 

hence we proved the inequality for j + 1. 

The inequality is true for j = 1, hence the end of the proof of Lemma [H where 
we may choose the value of R for A 6 [575] as R — x- Fhially we have the 
equalities, for all y such that > £o : 

U + (y,e) = 1 + (1 - e( |» (1 - A) E,>i 

= 1 + (1 "T) (1 " A) ((l(flF) ^>o ai,,(, c 0v))( T ^ F ) i 
V+(y,e) = A - 1 + (1 - A)(l - £(„)) 

= A - 1 + (1 - A)(l - S(tf))(^7)Ei>o W^X^?'- 
Using the estimates (fl8|) and the change of variable (fT6|) . for £ < R~ 1 the series 
S%(^x)C J is normally convergent and the following functions are well defined 

f ^(i/, g ) = i+ (i -y» ^(c, g ) 

\ y(y, e) = A - 1 + (1 - A)(l - £(y))CB(C, e)- 

It is straightforward to check that U and V solve system (fTU)) and that we 

have, for > £ , C(£) < |p hence foi@ e < e = |% and > Co we 

have U(y,e) = U+(y,e) and V(y,e) = V+(y,e). We extended the solution 
constructed for S [£o, 1[ to the region £ < -4. This proves Proposition [21 

4 Note that these inequalities depend on a given arbitrary £0 > 0. 
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2 The solution in the low density region 
2.1 Construction of the bounded solution 

In this section, we obtain the family of solutions of ([9]) bounded by \y\ A e ey when 
y — > — oo, that is in the low density region £ — > 0. Introduce the new variable 
f = — ey. Commonly, I call this solution the hypergeometric solution, because it 
has been observed that, in the model case po(%) — {—x — studied in [3] as 
well as in [B], the Rayleigh equation rewrites as the hypergeometric equation. 
Introduce 

r(s, e) = - A(£(-£))(£(-£))-i = £(i - £) = 1 + ei S '( S ,e). 
as e e e vs 

Wc define the operators R ei K e and through 

Re{9){s) =[ r r(y,s)e- 2 y(a-^))- X g(y,e)dy]e 2 %a-i))\ (20) 

J s e ° 

1 — A 2 /"+ 00 

K e (g)(t) = (1 - \)K*(g)(t) = — ^ y (21) 
These operators rewrite 



/■+oo 1 

fl £ (5)( S ,e) = / (-+e^'( 2 /,e)) e - 2 ^ s ) S -^^exp( £ ^A(5(y)-5( S ))) 5 ( 2/ , £ )dy. 

1 _ \2 p+oo i 

ife(ff)(*,e) = — : — / (-+e»S'(s,e))R e (g)(s,e)d S . 
We have the inequalities, for g uniformly bounded, (and A < v, which implies 

a-i r~ x < a~ir- x for t> s ) 

r+oo i p v 

\Re(9)(s)\<\\9\U -C-Hl-0e- 2y dy}e 2s e<\\9\\°o- (22) 
J s e e 

|\2_i| roo pv |A 2 — 11 f v 

\K e (g)(t)\ < !—-L Moo / T (s,e)^ds < L-—L\\g\\ O0 S-. (23) 
4 J t ewe 

Moreover, the following inequality is true: 

\g(s,s)\ <C p (^y=> \K e (g)(t,e)\ < ^-A- )Cp {^f + \ (24) 

In a similar way, we introduce 

K$(g)(t) = ^f t + °°±R>(g)(s)ds 
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Let £ > be fixed and < e < £q. Under suitable assumptions on g (we can 
for example consider g in C°°([to, +°°[) such that \d p g\ < C p y a ~ p for all p), the 
operators K e , R £ , Kq, Rq are well defined. Moreover, one proves that 

g(t,X,e) = J2^ n \m,e) (25) 



g (t,\) = Y,K { Q n) m) (26) 

are normally converging series on [to, +oo[, and that we have: 

g = 1 + K £ (g),g = 1 + K (g ). (27) 

Moreover, we know that g is defined on R, because the series ^2 — p \ ( £ ) p 

converges and is majorated by exp(|A 2 — 1|A^-), from the inequality (|24|) . We 
obtain the inequalities 

|. 9o (t,A)| <exp(^-=^),| 5 (t,A,e)| <exp(^l-^r 1 ). (28) 



We cannot thus consider the limit ( — ► in the equalities containing g as (|28|). 
We shall assume that A belongs to a compact set and that A > \. We prove 

Proposition 3 Let g be defined through \25)) . The family of solutions of the 
system Ul\) on (u,w) which is bounded by \y\ A e ey when y — > — oo is given by 

u(y, e) = C(F(t, A, e)+G(t, X, e)), t(y)wfa e) = v(y, e) = C£(y)[(X-l)F(t, A, e)+(\+l)G(t, A, e)] 

where C is a constant, t € [to, +oo[, t = ~ey and F and G are given by equalities 

CS3) and ([HP below. 

We have the estimates, for t e [to,£[ 

\g{t,\,e) ~ g (t,\)\ <C ei\go(t,X)\ 

\u(—,e) - u (--,e)| < C S"\uo(--,s)\ 

e e e 

\v{—,e) - v (—,e)\ < C £-|fo(--,e)| 

£ £ £ 

proof The system (| 1 1 1) rewrites on F and G given by Proposition [3j 

F'(t, A, £) = F(t, X, e) - |(i + eiS'(t, e))[(A - l)F(t, X, e) + (A + l)G(t, X, e)} 
G'{t, X, e) = -G(t, A, e) + + ehs% e))[(A - l)F(t, X, e) + (A + l)G(t, X, e)}. 

(29) 

A non exponentially growing solution of the system (|29|) is obtained through 
the following procedure. We denote by g(t,X,£) the function 

g(t,X,£) = G(t,A, £ ) e t (£(-f))^ i (g)-^ 

= G(t,X,£)e t t- 1 ^exp(~£i^S{t,£)). { ' 
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We first get, from the fact that F is bounded when t — ► +00, that 

F(t,X,s)e-H^e ei ^ ls ^'> = F(t, A, ^-^(-i))^ (f 

= A±i J t +00 (i + eiS'(s, e))s^e si ^ s ^e~ s G(s, A, e)ds 

(31) 

We deduce from the system ([25]) the equality 

j f (G(t, A, £ )e*r ^ exp(-e^ ±±±S(t, e))) = *zJl(l+ e Ss')e**-^ exp(-£^ e))F(t, A, e). 

Under the assumptions 17 bounded and satisfies the condition 

lim t _Kx,5(t, A, e) = 1 (32) 

one gets the equality 

g(t,X,e)-l = K e (g)(t,e). (33) 

Using the usual Volterra method and inequalities (f23|) , (I24| and (|28|) , we deduce 
that the only solution of (|33|) satisfying assumptions (|32|) is given through (|25|) . 
One gets G through (|30|) then F thanks to 

F(i, A, e) e - ( ^ (-)^ = (£)* ^- / r(s, e ) e - 2a r A s(s, A, e)ds. (34) 
v v 2 J t 

The first part of Proposition [3] is proven. 
Denote by (uq, Wq) the leading order term in e of (u, w) when t and A are fixed. 
Introduce Fo(t, A) and Go(t, A) through the equalities 

uo(t, A) = F (t, A) + G {t, A), A) = (A - l)F (t, A) + (A + l)G (t, A). 

The functions (Fo(t, A), Go(t, A)) are solution of 

f A) = F (t, A) - ^F (t, A) - ^G (t, A) 

I A) = -G (f, A) + ^F (t, A) + ^G (t, A). 

The second part of Proposition [3] comes from the following estimates on the 
operators R £ and K £ , valid for £ < £0 and t > to > 0: 

\Re(f) ~ i? A (/)l < C^\R x (f)\, \K e {g) - X A (.g)| < C 2 ei \K^(g)\, (35) 

from which we deduce the uniform estimates for g given by (|30|) solution of (|33| 

|fl(t,A,e)-flb(t,A)| <C 3 ei\g Q (t,\)\,t>t ,e<e (36) 

because the Volterra series associated with Kq is normally convergent in [to, +00 [. 
This ends the proof of Proposition [3] 
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Note that the previous estimates, as well as the behavior of the solution and the 
operator Ro, are valid only for to > 0, because, for example, Rq(1)(s) ~ 1 when 
s — > 0. The integral defining Kq is nevertheless convergent at +oo, because for 
t > to we have the equality 

/>oo -i 

2^(1)00 = 1 - / -{ V -)^e- 2 ^- s Uy. 
Js V s 



2.2 Construction of the hypergeometric solution for e = 

We prove in this Section 

Lemma 3 The solution (F (t, A), G (t, A)) constructed through fifflfy, fifflj), 

for e = is given by 

F (t,\)=e- t (U (t,\) + ^(t,Xj) 
Go(t,X)=e- t (U (t,X)-^(t,X)) 

where Uo(t) = 2~~i^~U{— 37-, — ^, 2t) the function U(a,b,T) being the Loga- 
rithmic Kummer's solution of the confluent hypergeometric equation (see \Tj). 

This allows to obtain the limit of the (Fo(t, A), Go(t, A)) for t — ► 0. 
The equation satisfied by Uo(t, A) = uo{t)e l is 

iltf - (2i + h)Ud + - 0. (37) 

Introducing T = 2t, we recognize (see [1]) the equation for hypergeometric 
confluent functions for b = — ~ and a = 

The family of solutions of this Kummer's equation is generated by two functions 
M (a, 6, T) and U(a, b, T). Note that T x ~ h M (1 + a - b, 2 - b, T) is also a solu- 
tion of (|3"T|). independant of M (a, b, T), hence U (a, b, T) can be expressed using 
M (a, b, T) and r 1_6 M (1 + a - b, 2 - 6, T). The family of solutions of (37]) which 
go to zero when T — > +oo is generated by U(a,b,T), called the logarithmic 
solution. It is the subdominant solution of the hypergeometric equation. 
The expression of the subdominant solution U(a,b,T) is the following: 

TTt hT )- - f M ^ T ) 1 _ 6 A/(l + a-b,2-fe,T) 

1 ' ' ' Bin7r6T(l + o-6)r(6) r(o)r(2-6) J 

where T is the usual Gamma function (T(s) = J °° t s e~ t dt). The relation be- 
tween J7 (a, b, 0) and (7'(a, b, 0) characterize the subdominant solution of the ordi- 
nary differential equation, and this particular solution has been chosen through 
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the lirniid when 



z — ► +00: 



U(a, 6, 0) = r ^ ( j Q h ] b) , \im z ^ +00 z a U(a, b, z) = 1. (38) 

As we imposed that g(t, e) — > 1 when t — > +00, we get that Go(i, A)e*i~^7r- _^ 1 
when i — * +00 and that there exists a constant C such that Fo(t, \)e t t 1 ~^ r — ► 
(7 when £ — * +00. Hence 

(F Q (t,\)+G Q (t,\)yt~^ -> 1. 

As T a U(a,b,T) -> 1, we get that r^rf7(-^,-i,2i) -> 2^. We thus 
obtain the equality 

r^U (t, A) = r^e^Fofa A) + G (t, A)) = 2-^r^U(-^^, --, 2t), 

2v v 

hence 

U {t,X) = 2- h ^U{- 1 -±±--,2t). (39) 
2v v 



Introduce 



We get that u (*) = 2"^C/(- i 2rr, -±, 2i)e"*. As w = Au - ^ = ((A 

dUg 

dt 



l)Uo — ' one deduces 



G (t,X) = ([/ (t,A)-i^(t,A))e- t ,F (t ) A) = (I7 (i, A) + ~-{t, AJJe"*. 

(41) 

Using [l] and (j37|) . we finally obtain 

G (t, A) -> 2-^(7o(A),Fo(t, A) -» 2~^G (A) when i -> 0. (42) 
We deduce the equality 

1 *+i r(i ~\~ ~) 1 A-i , 1 , A-i 

U o (0,X) = 2 X ~— — ^ hm^ +co t--- — [/ (t, A) = 2-+—. 

Note that we can deduce the expressions of F + Go and of Go- We thus 
check that 

(F + G )(t, A)e* = C (M(-i±*, -~ 2t) - C,(2t)* +1 M(1 + 2 + 1, 2t)) 

2 ^ 2^ v 

(43) 



5 r(i-t) E 

r(l + a-6) sin7r6r(6)r(l + a-b) 
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e*G (t, A) = C ((M - M') 2t) - C,(2t)*+i(M - M')(l + 2 + 

-CbC2^(l + ±)t»M(l + ^ 2 + i 2i)). 

(44) 

We note that (M-M')(-^,-±,0) = We deduce that e t (F +G o )(0, A) = 

2!-^C (A) and e*G (0, A) = 2-^C (A)(l-A), hence (A-l)e*(F +G )(0, A)+ 
2e*Go(0, A) = 0. In the next Section, we combine the results of Section [1] and 
of Section [2] 

3 Precise calculus of the Evans function. 

The Wronskian is related to a function independant of the variable t, called 
the Evans function, introduced below in (|4"B")) and denoted by Ev(X,e). In the 
present Section, we shall identify the leading order term in e of the Evans 
function, and all the terms of the form e^(A — 1) of the Evans function. We 

2 

shall finish by the calculation of the term of the form £>. More precisely, we 
prove 

Lemma 4 The function 

Ev(\,e)=£(y )W(yo) (45) 

is independant of y$. It is analytic in A and in si, e. Moreover, one has 
Ev(l,e) = 2(£)£ and d x Ev(l,0) = 2 1 ~^r(l + ±). This function is called the 
Evans function of the equation fP|). 

Using the expressions of ■^-(U{—ey,e)) and -£ju + , we have 

eW{y, e) = u + (y, s)(-sXU(-ey, e) + eW{-ey, e)) 
-U(-sy,e)(-eXu + {y,e) + -^v + (y,e)) 
= ^){£,{y)u+{y,e)W{-ey,e) - U(-ey,e)v + (y,e)). 

Hence we have the following constant function to study, which depends only on 
A, e: 

Ev(X,e)=£(y)W(y) = [Z{y)u+(v,e)V(-ey,e) - £(y)v+(v,e)U{-ey,e)]. 

We shall use the equalities, valid for all yo (and to) such that both solutions are 
defined (which means y € [--j^tr, -^g]) 

Ev(X,e) = e(» )W(»b) = 

We begin with the 

Lemma 5 27ie Evans function has an analytic expansion in X, which coeffi- 
cients depend analytically on e and £>. 
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For the precise study of the different terms of Ev(X, e), we introduce 

£ - £(-|), C = p Co = vt = C(t, 0), for i > i > 0. 

We check that the function Ev(X, e) is analytic in A and has an analytic expan- 
sion in e~s and e thanks to the equality 

^ 1 1 1-^-H _ 1 

which implies that the relation between i and £ is analytic in e and e», 
Assume from now on A > \ and v > 2 and replace £(y) by £>£". Using this 
Lemma, there exists two functions Bq(s) and Co (A, e) such that 

£7v(A,e) = £u(l,£)+£ (e)(A - l) + C (A,e)(A - l) 2 . (46) 



Direct relations Considering the limit in ([521 for e — 0, we obtain 

£u(A,0) = (A — l)e*(F + G )(t, A)[— 1 + (1 - X)vtA{vt, 0) + z/tB(W, 0)]. 

As this quantity is independant of t, we consider the limit when t — > 0, hence 
we deduce that 

Ew(A,0) = -(A- l)2 1 -^ i C (A). (47) 
Remark that this implies the identity 

(X-l)e t (F o +G )(t,X)[-l+(l-X)utA(vt,0)+vtB(vt,0)] = -2 1- ^Co(A)(A-l) 

(48) 

which rewrites 

e*(F + G )(t, A)[-l + (1 - X)vtA(vt, 0) + i/tfl(i/t,0)] = -2 1 - 2 ^C (A). 

In a similar way, we check that, for A = 1, U+ — 1 and V+ = 0, and g(t, 1, e) = 1, 
which implies 

G(t,l,e) = e- t (-)-ic 1 =e- t (^ (49) 

£ V 

from which one deduces 

Ev(l,e) = 2e t G(i,l,e)£ = 2(-)» . (50) 
From ([46]) , the unique root A(e) of Ev(X, e) in the neighborhood of A = 1 satisfies 

£to(l,e) 



X(s) - 1 



S ( £ ) + C (A(e),e)(A(e)-l)- 



The two first terms of the expansion of X(e) — 1 in terms of £> under the 
assumption v > 2 are thus given through 
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w s -,_ J5u(l,e) . 2 

[£ >~ ~~B Q (e) + Co(l,0)(\(e)-l) +0[£ " ) - 



As A(e) — 1 = — E g^ty + o(e » ) , we write 

. . _ .. E <l' e )... „ . + nf^ 

(51) 



A ( £ ) 1 ~~ Bofe)-Co(l,0)(Bo(0))- 1 Si;(l,e) + °( e " ) 



Bo 00 



0)(l,0)(Bo(0))- 3 (St;(l J e)) a + o(e5). 



One is thus left with the calculus of Co (1,0) and of Bq(s) up to the order 
1. For the computation of Bq(s), we need the behavior of the solutions of the 
overdense system for A = 1. 

As in SectionCQ we introduce o,(£) = a°+£a} + 0(£ 2 ) and bj(g) = ■ V>\ I 
O(e)- We recall that (A((,e) = E°li%(0C J and CB&e) = E°lAW- 
Introduce 

«(0 = E ?~ lb p v (0 = E e'^MO = E *(0 = E c j X- 

j>l 3>1 j>1 j>l 

Lemma 6 The following relations are true 

e\F + G)(t, 1, e) - e*(F + G )(t, 1) = + 0(e* ) 

2e'G(t, 1, e) - 2e*G (i, 1) = -2^ + 0(eS) 

C(t, e )_Co(t) = -^C 

CA(C, - m) - CoA(Co,0) = CC[fc(C) - w(() - ^((w'(() + w(0)} + 0(ei) 
CB(C,e)(l - ^y)) - Co^Co, 0) = eCb(0 - «(0 - ^t(C«'(0 + «(C))] + O(e^) 

For i/ie computation o/Co(l,0), one /ias 

e t (F + Go)(t,A)- e *(Fo + Go)(i,l) 



C (1,0) = -Zim, 



A->l,t->0 _ 



A — 1 

From these two results, one obtains the following 

Proposition 4 Introduce the function Ro(t) — ^- \nt — Kq (l)(t) — ^Bo(0)t~i , 
where Kq(1) has been introduced in &21\) and note that the terms B and Cq 
which have been introduced in \J^0^ are calculated through 



B (0) = -2 / sh- 2s ds = -2~^r(l + -) 
Jo v 



We have 



B (e) = B (0) + -^-j + 2(^-)i lim t ^R (t), 



and 

C (1,0) = J +o ° sherds -1 f+°° In se- 2s ds + ± /+ 00 ^-^-^(lXaJds 
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Let us begin with the proof of Proposition 2] We rewrite the Evans function as 

Ev(X, e) = [(A - l)e'(F + G)(t, A, e) + 2e'G(t, A, e)](£(j/) + (1 - A)(l - OC^(C, e)) 
-e'(F + G)(i, A, e)(A - 1 + (1 - A)(l - £(y))CB(C, *))■ 

(52) 

Remember that we have 

(A - l)[B (e) + G (A,e)(A - 1)] = Ev(X,e) - Ev{l,e). 

We thus deduce the equality 

5 (e) + C (A, £ )(A-l) = ^)^Mz^i£) 

-(1 - £(y))[e*(F + G)(i, A, e)(l - CB(C, e)) + 2e*G(t, A, e )CA(C, s)] 
+(1 - A)(l - f^JCAfCeJe^F + G)(t, X,e) 

Recall that G(t, A, e)e t = (§)"^9(t, A, e) and use <?(i, 1, e) = 1. We use also the 
relation (j34|) to get 



(F + ^G)(t, A, e)e-\^ = 2 ^(^ e -^g(s : A, £ )ds 

-^/rr A (f)^- 2 i|-2( 9 -i)]& 

(53) 

Note that we need two terms of G and of F + G, and that we use 

f = (1 - A )^A A (3)), 3 - 1 = (1 - A)K £ A ( 5 ). 
This will contribute to the term in C. Rewrite the first term of (1521) as 



* CM ,«)-*c(u.«) _ 2(£)i[ ©Z^i _ (C)^,,. 

A— 1 t/ A— 1 

Its limit when A goes to 1 is 2(~)* [-L ln(^) — K\ (1)]. Hence we get the identity 

B„(e) = -(1 - e-C-'JIe*^ + G)(t, 1, e)(l - CBi(C, e)) + 2e 4 G(i, 1, e)C4i(C, e)] 

(54) 

and the right hand side is independant on t. Using Lemma |51 we obtain 

B (e) =2^)i[^ln(§)-^(l)(i)] 

+ (1 - e ^C--)(-iHi)(l ~ CSi(C, 0) + 2CA 1 (C,0) + o(ei)) 

+ (1 - ^C'MFo + G )(t, 1)(1 - (B 1 (C, e)) + 2e*G (i, l)C^i(C,e)] 

from which one deduces 

B„(e) - 2(£)*[£ ln(§) - - ^(1 - CBi(C,0) + 2^,0) + (e*)) 

-(l-^C-')So(O) 

-[e\F + G )(t, l)(C(Bi(C. 0) - Bi(C, e)) + 2e*G (*, l)C(^i(C, e) - CMC, 0))]. 
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Using the relations G (t,l) = (%)"e~* and F (t,l) = 2e* J t + °° s^ e - 2s ds, one 
deduces that Go(t, l)and F (t, 1) goes to a constant when i — ► 0. 
Hence one gets 

B (e) = B (0) + (~)-lim t ^ [~ Ini - 2# 1 (l)(t) - Hb (0))]. (56) 

The second part consists in the calculus of Co(l, 0). 
Considering now e = in (|52p . one obtains the two identities 



B o (0) = -e\F + G)(t, 1, 0)(1 - Co£i(Co, 0)) - 2e 4 G(i, 1, 0)Co^i(Co, 0), Co = vt. 

B (0) + C (A, 0)(A - 1) = -e*(F + G )(t, A)(l - Co5(Co, 0)) - 2e'G (t, A)C A(Co, 0)) 

-(A-l)CoA(Co,0)e*(Fo + Go)(t,A). 

Hence 

C (A,0)(A - 1) = -(A - l)CoA(Co,0)e*(F + G )(t, A) 

+e<(F + G )(t, 1)(C5(C, 0) - CBi(C.O)) 

+ (1 - CoS(Co, 0))(e'(F o + G )(t, 1) - e 4 (F + G )(t, A)) 

+2e t G (i,l)(CA 1 (C,0)~C^(C,0)) 

+C A(Co, 0)(2 e 4 Go(i, 1) - 2e*Go(t, A)). 

We get (as we work for e = 0, we should write Co but we drop this notation and 
we use C = vt) 

G o (A,0)= -CA(C,0)e t (F a + G a )(t,\) 

-e^F + G )(M)C Bl(C l-f (C0) 

-(i - cb(C o)) et ( J '° +G °)( t ' A )- et ( Fo+G, °^ t ' 1 ) 

-2e t Go(t,A)C A(C,0) ^f (C ' 0)A " 1 
-(A(( o) 2etG o( t ' A )- 2e ' G o( t - 1 ) 

In this equality, one only needs the value for A — > 1, and it is independant of £■ 
We thus consider the limit when A — * 1 and C — * 0, hence one obtains 

r n n , v e t (Fo + Go)(t,A)-e t (F +Go)ft ; l) 
Go{l,0) = -hm c ^o,A^i — ■ 

Equality (|53|) rewrites 

(F+^G){t,\e)e-\^ 

Hence, considering the limit e — > 0, one obtains 

(fo + ^lG )(t,A)e- t (§)^ i 

+ (A"l)^r(§)^^o A (5o)e- 2 1d S 
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The value for A = 1 is thus (F + G )(i, l)e"* = 2 J^°° e~ 2s ds. Hence 
e-*[(^G )(t, A)(§)^ + (F + G )(f, A)(§)^ - (F + G )(t, 1)] 



1+A f< 
2A Jt 

Dividing by A — 1, one deduces 



+(^-i)^rc^[^(.9o) e - 2s Ks 



[-iG (t, A)(§)^ + (^+Go)(t,A)-( J b + Go)(t,l) ( ^ ) ^ + (Fq + Go)(t> 1} (i) 



A- 



Ml J- (fO^V^ 

+ i *rr(§)^^o A (3o)^ 2 1^ 

We consider the limit when A — > 1, and recalling that for e = one has j = s, 
denoting by JT(t, A) = ( f °+go)(M)-(Fo+G )(t,i) ^ we obtain 

e-*[-|G (t, 1) + H(t, 1) + (F + G )(t, 1)£ la §] 

2 r°°ln S e- 2s d S 
+ r«U^o 1 (l)e- Sto ]A 

Using again the integration by parts on the last term hence one gets 

e-*[-|G (t, 1) + H(t, 1) + (F + G )(i, 1)£ In §] 

2 /~ lnse"^ - ti^(l)(f)e- 2t - J, 00 £(si)Ki(l)e- 2 °ds 

We notice that the function R (t) = -± In t + 2K%(l)(t) + 1~ i B {0) has a finite 
limit when t goes to zero, according to ([56| . We have the equality K^{t) = 
\Ro(t) + i lnt - iS (0)ri. We deduce that 



// A (s i)^i ( i )e - 2sds = £ * (^)[Ii? ( S ) + ^ Ins - IS (0) S -i] e - 2 ^ S 

= vi- 1 // s i-i[Ii? ( s ) + -L ins _ lB o (0) s -^] e - 2 ^s 

In this last term, the only term which matters when t — !• is the term 

-iso(O)^- 1 J s- l e- is ds = -is (0)^-i[jT 6 2S ~ l ds - lnt]. 
Note that B (0) = -2 J + °° sie~ 2s ds. One obtains 



Bo(0) 



/"OG 1 

- / U»e- a da = -2-iT(l + -). 
Jo L v 
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3.1 Reduction of the Evans function 

Lower order terms Recall that the operator K^(l) is defined through ([2"T]) . 
We prove the following lemma of reduction: 



Proof of Lemma [6] It is enough to prove that the relation giving ( is 

1 1 

~vt? [y - 1)^ 

hence we deduce 



V V —1 

We thus obtain t = — , hence 

v v — 1 

We deduce that w(Q - w(( Q ) = (C - CoV(Co) + 0({( - Co) 2 ), hence w(() - 
w(Co) ~ (C ~ CoK(Co) = 0(f) and w(Q - w(( ) - (C - CoW(C) = 0(e#). 
We use e t G(t,X,s) — r s(t,A,e), hence for A = 1 we obtain 

e*G(t,l,e) = t? -1 . 
The equality giving F(t, 1, e) being 

/ + OG /" + 0O 

r(s, e)e~ 2s r]{s, e^ds = ~e- 2t T]{t, e)~ 1 +2 J e~ 2s f?(s, e) 

one obtains 

r°° 1 1 

e t (^ + G)(<,l,e)-e t ( J F 1 + G )(t,l) = 2e 2t / e~ 2s {- - - r - R r)d S . 



'T]{S,S) T)(S,0) J 



Similarily 



e t (2G(i,l,e)-2G (i,l)) = — -^[1 



77(5, e) 57(s,0) 
Using the relation 



ri(t,e)" y v-\ n ' ' " r](t,0) v 

one obtains 

^1-1 = ^^,0) + 0(0. 

»7(t,0) f-1 
This gives directly the two equalities of Lemma [6] 
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3.2 Limit for large k of the growth rate 

Recall that was proven in [TU] the following estimate on any value of 7 such that 
there exists a solution of © associated with A = ^| and e = kL : 



7 — * A = < /— w 7 -r — — when fc 



If we compare with (|13j) . one may see the difference between the result for 
Lq — > when fc is fixed and the result for Lq > fixed and fc — > +00. Note for 
example that the limit of 



gk 



^i + (M*)*r(i + i) 

when /c — > +00 is +00 because v > 1. This is not surprising because we did not 
get the lower order terms up to the order e of the expansion of A. Remark that 
the term in e comes from the terms in e in the functions A{C,,e) and B(£,e). 
We have the following result (according to [TO] ) 

Lemma 7 There exists fc* > such that, for all k > k*, there exists a realj(k) 
and a non zero solution u(x)e lky+1 ^ t of the Rayleigh equation such that 

- < 7 (fc) < A. 
We have the following behavior of the eigenmode 

\\Po u W + \\Po u '\\ + IMI + + IK'll < +°° 

As the result of this Lemma is important for the nonlinear analysis, we rewrite 
an idea of the proof, based on Remark 8.1 of [10]. We denote by L 2 x the space 

Po 

1 

of functions u such that p^u € L 2 (R) . 

Finding 7 is equivalent to finding as an eigenvalue (in L 2 (R)) of 
1 -3 d i d ~h . 1 9 . f x 

~Y 2 ^ Tx {p0 Tx p °~ ) + l ~7 M ) - 

This operator rewrites — + 1 — -^k (x) + k~ 2 Wo{x) where Wo (a;) = 

ifcg(x) + |(fco(x)) 2 , which is bounded when^- is bounded (or equivalently when 

fcg is bounded). We introduce the operator Q = — p-p 2 j^.{poj^Pq 2 ) + which 
is coercive, thanks to the Poincare estimates, for fc large enough. The eigenvalue 
problem rewrites 

2 

— G a p (Q-h Q-^). 
9 

Under the (natural) hypothesis that fco has a nondegenerate minimum Lq, one 
deduces that for k large enough one has at least a value of j(k) such that 
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Lq < ( 7 (|))^ < 4io using usual results on semiclassical Schrodinger operators 
which potential has a well. 

We thus constructed v € L 2 (R) and j(k) such that v is the eigenvector of 
Q~^koQ~^ associated with the eigenvalue ilML. 

To v is associated a solution of ([4]) which is it = p 2 Q _5 f, it' € £ 2 1 , u 6 L 2 1 
Remembering that it solves 



Po PS 



-u" + k 2 u - k (x)u' - y4rW u = °J 

multiplying this equation by it and integrating, one gets 

/— — 9k 2 — — 
k {x)p 2 u.[j-y-^p§u + p§u]dx 

hence, using the hypothesis 

P0P0 1 < M 

one obtains (the norm on the Sobolev space H 1 is \\u\\\ = j(u') 2 + k 2 u 2 dx) 
\\u\\i<M[j^— 2Mu\\ + \\p>u'\\} 

1 

hence a control on the H 1 norm of u (instead of having the weight p§ ). 
Moreover, as it" = r-rcyp ko(x)u + k^(x)u' — k 2 u, one deduces that it" S £ 2 , and 
we have iteratively the control of it in £P (s < s ma x, according to the number 
of derivatives of fco that we consider). 
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4 Towards a non linear analysis 



We show in this Section that the result of Guo and Hwang [8] can be extended 
in our set-up, even if the density profile po(x) does not satisfy the coercivity 
assumption (3) of [5]. The quantity k (x) — plays a crucial role. It has a 
physical interpretation, being the inverse of a length: it is called the inverse of 
the density gradient scalelength. We need the assumptions 

_ i 

(H) ko(x) bounded ,ko(x)p 2 bounded. 

Note that ko bounded is fulfilled in the case studied by Guo and Hwang (where 
po is bounded below), and in the case of the striation model (studied by R. 
Poncet [H]) but is not automatically fulfilled by a profile such that po(x) — > 
when x — > — oo. However, for the particular case of the ablation front profile, 
we have k (x) = Lq £(j^) v (l — £(7^)), hence it is bounded and belongs to 

Before starting the proof of Theorem [2J which is rather technical, let us de- 
scribe our procedure. 

Firstly, we prove that the linear system reduces to an elliptic equation on the 
pressure, from which we obtain a general solution. We identify a normal mode 
solution of this system using the first part of the paper. 

Once this normal mode solution U is constructed, with suitable assumptions on 
the growth rate, one introduces a perturbation solution of the nonlinear system, 
which initial condition is SU\t=o and an approximate solution V of the non 
linear system which admits an expansion in 5 N up to the order iV with the same 
initial condition. 

Using the Duhamel principle for the construction of the j— th term of the ex- 
pansion in 5 of V N , one obtains a control of all the terms of V ■ 

The natural energy inequalities are on the quantities p§u J , p§v J , p 2 p J , p 2 p 3 . 
We verify that the properties of po(x) imply that we can deduce inequalities 
on v? , , Pq 1 p> and T- 7 . 

Note that we have, as a consequence of the method that we chose, a control in 
t s e At of the H s norm of all solutions of the homogeneous linear system (with 
any initial condition U(x, y, 0)), and a control by e J7 ' fe '* (with no additional 
power in t) of the H s norm of the j— th term of the expansion. 

Remark 1 When an initial value mixes eigenmodes, the H" norm of the solu- 
tion behaves as t s e At . If one starts from a pure eigenmode with ^ < j(k) < A 
the exponential behavior comes at most from the growth of the pure eigenmode. 
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4.1 Obtention of a solution of the linear system 

Consider the system 

d t <J + Po«i = fo 
Pod t vi + d x p = ag + f x 
Pod t v 2 + d y p = f 2 
d x vi + d y v 2 = 

i _ i 

We know that the relevant quantities are Pg^i^, Po 2 c, and we denote these 
three quantities by X, Y, t. To have the same behavior when po — ► 0, consider 
ip such that, once ip is obtained, we revert to v\ and v 2 using v\ = —d y (p 2 ip), 

_ X 

v 2 = d x (p 2 ip). Introduce 



b = Pa 2 [dy(povi) - d x (p v 2 )}. 
The system on vi,v 2 ,a,p implies the two equations 

f d t b = gd y T + 2 (dyfi - d x f 2 ) 
\ d t T + k Q (x)X = Po*f . 

We obtain ip from b through the elliptic equation 



We then revert to X through the equality X = 
is obtained through the elliptic equation 



-b. 



(57) 



(58) 



(59) 



-dyip. Finally, the pressure p 



Pod x (p 1 d x p) + d^p = p§ [p§ d x (p 2 r)g + p 2 d x (p 1 fx) + p 2 f 2 ] 
which rewrites 

i 1 i 

Ap - k d x p = p 2 [d x rg - -k rg) + p 2 (div/ - fco/i)] 

Hence we solve the system 



d t T = k d y ip(b) + p 2 f 

_ i 

d t b = gd y T + p Q 2 (dyfi - d x f 2 ) 
T (°) = T o(x, y), b(0) = b (x, y) 
- {\k' Q + = -b 

, Ap - k d x p = p 2 [d x rg - \k rg) 



(60) 



(61) 



Po 2 ( div / - Wi)] 



which has the same properties as the system (13) of [5], the Poincare estimate 
being still valid. 

From b and r, one reverts to X and Y, hence a solution of the system. Moreover, 
one checks that (X, Y) 6 L 2 (R) (according to the energy equality), hence X g 
H 1 (R) under the assumption fco bounded. 
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Proposition 5 Under the hypotheseses (H), and under the hypothesis hj € 
L 2 , j = 0,1, 2, the functions U\ , V\ , T% , p\ solution of 

d t Ti - k Q ui = h 
p d t ui + d x pi + pogTi = hi 
p d t vi + dypi = h 2 
d x ui + dyVi — 

satisfies ui(t),vi(t),Ti(t) G L 2 when it is true for t = 0. Moreover, one has 
Po V(i)eL 2 (R 2 ). 



Proof The proof of this result follows two steps: first of all the assumption fco 
bounded implies that p§ u\ , p§ V\ , p 2 Vpi , p§ T\ belong to L 2 . We thus multiply 
the equality dtU\ + /5q 1 Vpi + Tig = h by V(pq 1 p). We get, integrating in x, y: 

J (V qi ) 2 + k (x)q 1 Vq 1 .e 1 + T ig y qi = J hVqx 

from which one deduces 

||V<zi|| < max(fc Po H\phx\\ + <?l|7ilU + \\h\\- 

i _i 

It is then enough to use the Poincare estimate between p§ q\ and p 2 Vpi to 
obtain the estimate on Vgi , from which one deduces the estimate on q± . 
Finally, from the estimate on q\ and on V?i, multiplying the equation on the 
velocity by u\ and integrating, we get the Gronwall type inequality 

j t \\ui\\<C\\ qi \\^ + \\h\\ + ff ||T 1 || 00 

hence a control on ||ui|| on [0, T] for all t as soon as it is true for t = 0. 
The system writes 

d t T + u.S7T = uTko(x) 

d t u + (u.V)u + TVQ + TQk (x)ei = (1-T)g (62) 
divu = 

In the system (|62p. appear only quadratic terms. When one wants to deduce 
the term of order N in the system, plugging in the expansions T N , u , v and 
Q N one obtains source terms of the form 

Sn = Y^=2 UjT N _jk (x) - Ujd x T N -j - vjdyTN-j 

Ri,n = " Y^=2 Ujd x u N ^j + v.jd y u N ^j + Tjd x Q N ^j + TjQx-jk^x) 

R2.N = - Ujd x VN-j + v-jdyVN^j + TjdyQN-j 

and the system rewrites 

d t T N - u N k (x) = S N 
d t u N + d x Q N + Q N k (x) + gT N = R X>N 
dtVN + 9 v Qn = R2,N 
d x u N + d y v N = 0. 



(63) 
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Higher order Sobolev regularity (preparatory equality) One of the 

main tools that we have to use is the divergence free condition, in order to get 
rid of the pressure p or of the reduced pressure Q when obtaining the energy 
inequality. Recall that the system (|6"5)l rewrites 

d t T N - u N k {x) = S N 
p d t UN + V(poQw) + gpoT N ei = p R N 
divu/v = 

where R N = (i?i,jv, R.2,n)- 

Denote by Gn — PodtRN — ffPoSWei. Applying the operator d t d™ n to equation 
on the velocity and using the equation on the specific volume, one obtains 

d^( Po df 2 u N ) + Wd t d^( Po Q N ) + gd^(p k u N ) = d^(G N ). (64) 

One deduces the 

Lemma 8 For all n, one has the estimate 

MdldZ n u N \\ < C n (J2 \\pld p xP u N \\ + \\p$B£,G N \\). 

Moreover, as the coefficients of the system depend only on x, this inequality is 
also true with the same constants when <9™„ is replaced by d™ n dy q for all q > 0. 

Proof One notices that ([64]) writes 



n— 1 n— 1 



7 b — -L 7 b — 1 

p d? 2 d^u N +V(d t d^(poQ))+gko(x)p Q d^u N = Gn-J2 ClP^dl.dluN-gY, C nP 

Multiplying by d^ 2 d" n UN and integrating, using the recurrence hypothesis that 
\\pld^d p xP u N \\<C p ( l|pj^™^||)+2g 2 A 2 ( ]T \\pl^u N \\)+\\G^-i 

m<p—l m<p—l 

as well as the inequalities 

\ko(x)g\ <A 2 ,| Po -V^| <A P 

(which are true as soon as fco is a C°° function which derivatives are bounded, 
because p f — fcoPo ) one obtains the inequality 

\\ P ld^u N \\ < c n (J2 WpI^™ 



Lemma [8] is proven. 
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4.2 The energy equalities 

Note that the system for the leading term of the perturbation is the system (p5t 
with a null source term. Owing to this remark, we shall treat the general case 
and apply the equality to the particular cases. 

Multiplying by dtd" n UN and integrating, using the divergence free relation, 
one obtains 

/ d" n (pod^ 2 UN)-d^ n d t UNdxdy + J gd% n (pokoUN + poS^)e\.dtd^ n UNdxdy 

J d% n (podt S N ).dtB£n u N dxdy. 

In this equality, we can consider (for Sobolev inequalities) the term containing 
the largest number of derivatives of upf. We obtain, denoting by 

&n = d^(p d 2 2u N ) - po<9™ re <9 t 2 2Ujv 
B n = d" n (p k Q u N ) - p Q k Q d" n u N 

the equality 

/ pod% n d% 2 UN-d% n UNdxdy + J gp k d™ n UN-d t d™ n UNdxdy 
+ j .d™ n d t u N dxdy + J gB^ .d t d™ n u N dxdy 

J d™ n {pod t RN)-d t d™ n u N dxdy - / gd™ n (p S N )ei.d t d" n u N dxdy. 

The terms R% and contain only derivatives of order less than n — 1 , hence 
it will appear as a source term in the application of the Duhamel principle later 
on. The two first terms of the previous equality are the exact derivative in time 
of 

EnW = k[ Po(d r ^d t u N ) 2 dxdy+ f gp k (d^u N ) 2 dxdy}. 



2 . 

The energy equality is thus 

E%(t) = EZ(p)+ [ t g^(s)ds 
Jo 

where 

9n{t) = I d^ n (pod t RN)-dtd" n u N dxdy - J gd^ n (p S N )e 1 .d t d^UNdxdy 
= -(/ R% .d^dtUNdxdy + J gB% .d t d^u N dxdy). 

Note that this source term satisfies 

\g%(t)\<\\pld» n d t u N \\ L2 K»(t) (65) 

where one has 

K»(t) < \\Po h R%\\ + \\p^d-4 Po d t § N )\\ (6g) 

+\g\l\\Po"B%\\ + \\dUposm- 

We are ready to prove the Duhamel inequality associated with this problem, 
using gko(x) < A 2 . 
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4.3 The Duhamel principle 

Two versions of the behavior of the semi group will be deduced. The first one 
corresponds to the general case for the terms in S 2 at least. 
We consider the (general) system 

p (x)d 2 2 w + V(po<9 t (2) + gk p wei = M , divw — 0. (67) 

with the initial conditions 

^|t=o = 0,^| t =o = 0. (68) 

Note that this system is easily deduced from the system obtained for the N th 
term of the expansion in S of the solution. 

Proposition 6 Assume that there exists two constants K and L, with L > A, 
sich that 

\\p-~ 2 M\\<Ke Lt . (69) 

The unique solution of the linear system |d>7| ) with initial Cauchy conditions 
168\) satisfies the estimate 

< 2K ., c Lt 

\\pld t w\\ <^(1+7e4f) 4 ^ 
I Md%w\\< K(l + I ^ F )e Lt 

Proof We begin by multyplying the equation (|67|) by d 2 2 w and integrate in 
space. One deduces that 

\\pld 2 2 w\\<A 2 \\plw\\+Ke Lt . 
We will make use of this equality later. 

Let us multiply the equation (f6"T)) by dtw. We obtain the identity 

li\ J Po( d tw) 2 dxdy + i J k p k w 2 dxdy = J M(x,y,t)d t wdxdy. 

Integrating in time and using the initial condition (|68p as well as the estimate 
(|69p . we obtain the inequality 

p (d t w) 2 dxdy < A 2 ^ J k aPa w 2 dxdy + 2K J e Ls ||p| d t w\\(s)ds. 

Let us introduce now u(t) = J* | \pfi dtw\\(s)ds. We obtain, considering J p^uPdxdy, 
that ^ 

Mw\\{t)< [ \\pld t w\\(s)ds 
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that is 

Mw\\(t)<u(t). 

Hence the inequality 

{u'{t) f < A 2 (u(t)) 2 + f 2Ke Ls u\s)ds. 
Jo 

From this inequality, we deduce that 

(u'(t)) 2 < A 2 (u(t)) 2 + 2Ke Lt u(t) 

hence 

u'{t) < Au{t) + ^2Ke Lt u(t). 
Introduce h such that u(t) = (h(t)) 2 e At . We obtain the inequality 

2hh'e M < sFXKhe^ 1 



hence 
that is 

which leads to 



, , , 1 / L — A j 

ft' ft) < -VlKe—* 
v j - 2 



, , „ \J2K i-A t 

2K Tt 



The estimate on u'it) follows, using (v!) < A 2 u~ + 2Ke u. We thus, by 
integration, deduce another estimate on u. The estimate on pod 2 2 w is the con- 
sequence of (jSJ. 

If one wants a general formulation of the Duhamel principle (taking into ac- 
count non zero initial values), one states the following proposition, which will 
lead to the result of proposition [51 hence allowing a mixing of modes and a weak 
nonlinear result. The mixing of modes is not our purpose here, but we shall 
not speak of weak nonlinear results. See Cherfils, Gamier, Holstein [4] for more 
details. 

Proposition 7 The solution of 

1 d I [/„ IF, X \? PO „ U. 



, ,;( / iPo{dtu N ) - g—p {u N ) )dxdy) = g(t,x,d t u N ) 

2 at J p 

with initial condition 9 t itjv(0), ujv(O), with the assumption 

i 

\g(t,x,d t u N )\ < K(t)\\p§d t u n \\ L 2 
where K is a positive increasing function for t > satisfies the inequalities 

WpIunW^ < [d + J* V K(s)e-^d s}e^ 
\\pld t u N \\ < [Cx + /* ^K(s)e-^ds] 2 e At 
where C\ depends on the initial data. 
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Proof Wc deduce from the energy equality the following inequality: 

pa(x)(d t u N ) 2 dxdy-g / k (x)p Q (x)u 2 N dxdy < C ,++2 / K(s)\\p§d t UN\\L 2 { s )d& 



1 



where Co = / po{x)(d t UN) 2 (0)dxdy — g J k o (x)po(x)u 2 v (0)dxdy and Co,+ = 

i t i 

max(Co,0). Consider now the function u(t) = \ \p§ un(0)\ |+/ ||Po dtUN(s)\\ds = 

— t — - 

||PqUjv(0)|| + L \\p§dtUN\\{s)ds. We notice that u'(t) = \\p§ dtUN\\(t) hence 

u'(t) > 0. Recall that gko(x) < A 2 . The inequality implies 
{u'{t)f < A 2 (u(t)) 2 + C ,+ + 2 /„* K(s)u'{s)ds < A 2 (u(t)) 2 + C ,+ + 2K(t)u(t) 

s- IK i K(t)\2 . n K(t) 2 

< (Au + -±-±y + c .+ - -y-. 

Use now the inequality (a 2 + b 2 + c 2 )s < a + 6 + c for positive numbers a, 6, c 
to obtain 

u'(t) < Au{t) + + y/2K{t)u(t). 

Introducing v(t) = u(t)e~ At which satisfies v(t) > u(0)e~ At , we deduce 



v'(t) < y/Q^e~ M + ^f2K(t)e- At v(t). 
• Assume w(0) > 0. We obtain, denoting by h(t) = y/v{t) 



2hti < ^/C^e- At + ^2K(t)e- At h(t) 

hence 

2h'<{^)^e-^ + ^2K{t)e- A K 

We deduce the inequality 

h(t) < h(0) + A-i(^±)h {1 - e -¥) + _L J* ^K{ S )e- A s d s. 
which imply that there exists A and B such that 

u(t) < (A 2 e At + B2 e M (J o \jK{s)e- As ds) 2 ). 
• Assume u(0) = u'(0) = 0. As Co,+ = 0, we have the inequality 



u'(t) < Au(t) + v / 2K(t)u(t) 
from which one deduces, with the same notations as above, that 



h'{t)<\l l -K{t)e-^ 



32 



hence with h(Q) = one obtains 



h(t) < ^K(s)e- As ds. 
Assume finally u(0) =0 and u'(Q) > 0. We obtain 



v'(t) < TO^e- At + yj2K(t)e-**v(t). 



Introduce v(t) = v(t) — \JCq. + ^—\ — . We have 



v'{t) < ^2K(t)e-"(v(t) + ^C^LjLir) < y 2K(t)e-^(v(t) + 
from which one deduces the inequality 

J ^2K(s)e- As ds. 



A 



In all the previous cases, we deduced the inequality u(t) < [Ci+L \J K(s)e~ As ds] 2 e Ji 
Using finally the relation 

d i i 

-J t \\p0U N \\ < ||p 2 9 t WAr|| = u'(t) 

we get 

WpIUnW < u(t)-u(0). 
These are the two estimates of Proposition 

Of course, the proof is much simpler in the case we are interested in, that is 
dtUN — 0, un — 0, where (using the notations of this paragraph, Cq = Co,+ = 
it(0) = u'(0) = 0), where one deduces easily 



^ju{t)e- M < / 2~^K(s)e- As ds. 



4.4 H s estimates for a general solution of the linearized 
system 

The H s inequalities for the solution of the homogeneous system We 

consider the system satisfied by the leading order term of the perturbation of 
the Euler system (which is the system particular case of for N = 1. 
We prove in this section the analogous of the Proposition 1 of [5] , with a slightly 
better estimate which shows essentially that the relevant growth rate is, up to 
polynomial terms, A: 
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Proposition 8 Let Ti(t),ui(t) be the solution of the modified linearized Euler 
system \62\l . There exists a constant C s depending only on the characteristics 
of the system, that is of ko and g, such that 

ll4T l (t)||jr. + ||pfri(t)||jJ. <C.(l+t)'exp(At)(||p o i T 1 (0)||ir. + ||pJiZi(0)||ir.)- 

Note that in these inequalities (which are general) a power of t appears in the 
bound for the norm H s . This is the general case. Note that similar estimates 
were obtained independantly by R. Poncet [15] . 

An important feature of this result takes in consideration an initial condition 
which is not an eigenmode of the Rayleigh equation, and which is a combination 
of different eigenmodes. As we shall see in what follows, the interaction of these 
different eigenmodes lead to a linear growth of the form (1 + t) s e At for the H s 
norm of the solution. 

Proof We prove in a first stage the H s inequality result for the system satisfied 
by (Ti,ui,vi,Qi). We use the pressure p\ in the analysis. The system imply 
the equation 

p {x)df 2 ui + VdtPi = p gk' u\. 

We apply the operator D m p to this equation. The energy inequality deduced 
from 1|65|) and from the inequality ()66|) is 

{{uiyfKk^ulf+C + Kl^uiit) 

where we have the estimate 

Kl{t)<\\to*Q\\ + \9\\\p**Bi\\. 

1. principal term 

i 

The inequation on HpqWiII writes 

(||l4"ill) 2 <A 2 ll4«ill 2 + Co 
hence one obtains the inequality 



HaKII < |l4"i(°)ll cosnA *+ Y + ||p<Ni(0)ll 2 sinhAt < D e At . 

2. derivative of the principal term 

In the inequality obtained for D\. p ui, the source term g\ is bounded by 
M Doe At \\Di. p dtui\ \ because it contains only derivatives of order n— 1 = 0. 
We have thus the inequality 

{{u\)') 2 < K 2 {u\) 2 + Co + 2MD e At ul(t) 
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from which one deduces 



MD _ An2 „ .MD , 

,2„2At 



hence 



that is 



((u\y W) 2 < (A«J + ^p e A *) 2 + Co C-^ fe 



(uly(t)<Au\(t) + ^e At + VCo~ 
d ,-l„-AU / M£ >o , fpr At 



from which one deduces 

< A- 1 ( v ^+ML>oi + Au}(0))e At . 

3. Greater order term: 

We prove thus by recurrence that there exists A n and B n such that 

u 1 n (t)<(A n + B n t) n e At , 
according to the inequality 

^ foAMv-Ats s fFT -At . (A„_i + tBn-lf* 1 

5 («n(*)c ) < VC n e + ^ ■ 

One deduces the same inequality for 4?u\{t). 

4. In the derivative D np , the only term which matters for the order of the 
power of t is n, hence one deduces that 

J2 (WpI D„, p un\\ + \\pI D n , p d t u N \\) < (C s +tD s ) s e At 

n+p— s 

Proposition [5] is proven. Note that this improvement does not change the be- 
havior of the approximate solution we intend to construct, because for a normal 
mode solution 

u(x,y,t) = u{x)e iky+ ^ k) \ 

where r y(k) has been calculated and where u(x) is solution of the Rayleigh 
equation, one has the following equalities: 

\\plD m , P ui(t)\\ = MD m , P ui^)\\e l{k)t 

WDrn^m = \\T 1 (0)\\e'rW t (70) 

||4^m, P QiWH = ||pJi? mrf> Qi(0)||eTW*. 

Remark that, according to Lemma[7j and to the equality ikQi (x, y, t) = ^^-d x ui(x,y,t), 
we have also the relations 

\\D m , P Mt)\\ = ll^ m , P «i(0)||e^ t , . 

l|An*Qi(t)ll = IIA., P Qi(o)||e^ fe )'. 
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4.5 The H s inequalities for the linearized system 

We consider the system (|6"3"|) . We apply the operator D m ^ p = d™d?. This system 
becomes 

d t D m . p T N - D mtP u N k (x) = D m , p S^ + J2q=o C p D m, q u N k^ p) (x) 

dtD m ^ p VN + d y D m ^ p QiM = D m ^ p S^ 
v N = 0. 

(72) 

We notice that this system writes as the system (|63p with a source term involving 
derivatives of the solution at a lesser order of derivatives in x. 
We introduce 

<{t) = \\ P ld^u N (0)\\ L 2 + [ \\pld t d^u N ( S )\\ L 2ds 

Jo 

and v»(t) = u%(t)e- At . 

We are now ready to study the behavior of the lower order terms of the 
expansion, assuming that we found a j(k) such that ^ < j(k) < A. 
We have to deal in a second part with terms of the form u%, where N > 2. In 
this set-up one has to use Proposition [71 because we cannot obtain the sharpest 
inequality using the estimate v! < Au + + s/Cq. 



4.6 Inequalities for the following terms of the expansion 

Recall that from Lemma [Jj (proven in [10]), there exists a normal mode solution 
of the linearized system of the form u(x, fc)e lfca+7 W* where ^ < 7(/c) < A. 
With this normal mode solution one constructs an approximate solution of the 
nonlinear system, of the form 

T N (x,y,t) =l+£f=i<^ (x,y,t) 
u N (x,y,t) = £. j 6 j Uj(x,y,t) 

v N (x,y,t) = Y,j=\ 53v i( x ,yi t ) 
Q N (x,y,t) = Qo(x) + j:^ 1 ^Q j (x,y,t). 

There is an important Lemma, which depends on Hypothesis (H): 
Lemma 9 The functions Uj , Vj , Qj , Tj belong to L 2 . 

The proof of this Lemma is a consequence of Proposition O which will lead to 
the control of the source term of the linear system on Tn, un, vn, Qn- 
We shall use the estimates of Cordier, Grenier and Guo [5], and the method of 
Guo and Hwang [8 to give an H s estimate of T N , u , v N , Q N and a L 2 estimate 
of T N — Tq — STx, u N — Sui,v N — 6vi to obtain a lower bound on T N , u N , v N . 
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We prove in this section the H s estimate un in the weighted norm ||pq .||. Using 

_ i 

the assumption k a p Q 2 bounded, we deduce estimates in H s for u^. The first 
result reads as 

Proposition 9 There exists constants Cq and A p , depending only on the char- 
acteristics of the system (namely g, k$(x) and its derivatives) and on the H p 
norm of the initial data such that 

<(*) < {C p a ) N {A p ) N ^e Nl[k)t . 

Remark 2 This estimate relies heavily, as in JW, on the quadratic structure of 
the nonlinearity, and that we give the precise estimate on the constant Cj which 
appears in (13) of jTjj. This estimate could not be obtained in the set-up of Guo 
and Hwang /<§/ because the nonlinearity was written using puS/u, hence a cubic 
nonlinearity. 

A second comment is the following: the inequality ^(k) > A allows us to forget 
the coefficient (1 + t) s in the H s estimate for a general solution of the linear 
system (obtained in Proposition^ . This is a consequence, as we shall see below, 
of the relation 



At , 

e ; e 



1 1 



(N 7 (k)-A)s ds < 1 N-y(k)t 

~ Nj(k) - A 



(to be compared with the relation e At L e^ A A ^ s ds < te At ). 
Case N = 2 

Recall that we have the following system 

d t T 2 - k Q (x)u 2 = —Uid x T\ - vxdyTi - T±Ui 

Po(x)d t u 2 + V(/9 Q2) + gT 2 = -po{x)[ui.VUi] - poTid x Qi - p'qQiTx 
divu 2 = 

We have thus the estimates 

\\pld t S^\\ + \\ P ls]\\<C]e 2 '<^ t . 

This means that K$(t) < D 2 e 2 " i{k) \ hence 

\\pld t u 2 \\ < (Co + ! y/2fhe^-i>ds) 2 e M 
Jo 

hence the inequality 

ll4^ 2 || + ||Pckl| <M e 27Wt . 
We need to derive estimates for the terms T 2 and Q 2 . For the term T 2 , one has 

~dl\ [po T idxdy= J k (x)p u 2 T 2 + J S 2 p T 2 dxdy 
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from which one deduces the inequality 

J t \\pH\\ <M\\plu 2 \\ + \\p$Sl\\ < (MMo + ^V 7 ^ 
hence the estimate 

llp^WII < WpIt2(0)\\ + Cl ^ k) M0 (^ {k)t !)• 
As for the estimate on Q 2 , one deduces 

8x(fio 1 d x (p Q 2 )) + d 2 y2 Q 2 + gd x T 2 = divS 2 

which imply estimates on Q 2 - 
Case N > 3. 

We start with the induction hypothesis that, for j < N — 1, there exists Cq and 
A such that 

I €T, J 1 -l-l I ^ £^€7^ 1 1 -l-l I 1 1 -I- 1 ««3> 11 -I- 1 1 i»# 11 -l-l I /^^3f 1 1 < ^4. J "— x Cg e^-r 

and that the derivative in time of all quantities is bounded by jj(k)A : >~ 1 C 'e : ' 1< ' k ' >t . 
Thus there exists M (independant on the number of terms which appear in the 
source term and which depends only on the coefficients of the system) such that 
the source term of (|66[) for n = is bounded by: 

K^(t) < MA N - 2 C^N 2 - f {k)e N ^ t . (73) 

Note that in this estimate the N 2 term comes, one from the number of the terms 
in the expansion J^Lo* ^j^N-j and a second onc0 from the derivative in time 
which appears in the source term dtS N . We thus obtain, using 

h N {t) < J ^2K^(s)e- As ds 

the inequality 



h N (t) < J2MA N ~ 2 C^N 2 j(k) / y/eWiW-^'ds 

Jo 

which yields 

h (t) 2 < A N-i C N c (N^k)-A)t 8MiY 2 7 (fc) 
ftnW °° 6 (N^-AfA- 

The choice of A is thus induced by ^j^J^a — ^ ^ or a ^ ^ (forgetting that 
we have to be more precise to obtain estimates not only on ujv but also on Tjv) 



6 Note also that if we consider a cubic model, the number of terms in the source term is 
N(N — 1), hence adding a derivative in time we get N 3 in the estimate. As we can see in the 
following lines, this gives a less efficient estimate 
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hence the simplest choice is A = 8 f[?^_ k A ■ The value of Cn is thus given by the 

norm of the leading term (Ti, ui,«i, Qi). 
The final estimate is 

\\Po u n\\ < O A e w 
We proved the assumption (|T3")) . 

We use this result and the estimates for a normal mode solution (on which no 
powers of t appear for the norms of the derivatives). We obtain 

i Cn _i_ 1 I l /* JV-y(fc) — A 

"(0) Jo 

hence as Nj(k) > A one gets 

We deduce the inequality (using (a + b) 2 < 2(a 2 4- b 2 )) 

i 

u(t) < 2(M0) + u(0)) + A-^^nV* + 2( jV7( ^_ A ) 2 C jv7 (fc)eW. 

Remark If the system has a cubic source term, at each stage of the construc- 
tion one gets A1 N ~ 1 C N as estimate, hence the convergence of the infinite 
series is not ensured by these estimates. 

4.7 Estimates for the approximate solution 

In this paragraph, we derive estimates on the global approximate solution. We 
shall use throughout what follows the Moser estimates, that we recall here 

H^C/^IU^Cdl/IUIblU + llplUII/ll,) (74) 

and 

\\D a (fg) - fD a g\\ L i < CdlD/IUIIflll.-! + ||fl||oo||/||.) (75) 
and the Sobolev embedding H/H^ < C||/|| s for s > f and HV/IU < C||/|| s for 
s > | + More precisely, we prove that 

Proposition 10 For all 9 < 1 and for all t < In s J oA , we have 
n 1 „N u -~ SACoerW 



||T iv -1\\ H . + W^Wh* + ||Q iv - qo\\ B . < C 
\\T N - l|| La > WT^O^Se^ - AC*C 3 5 



1 - SACoefW 

e 7(fc)« 
1 - SACoe-tW* 



\\u N \\ L 2 > WuMW^Se^ 1 - AC 2 C S 6 



2 



e 



7 (fc)t 



1 - SACoeiW 
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p7(*)t 

\v N \\l* > \\vi(0)\\LiSe^ k]t - AC*C 3 5 2 



1 - <L4C eTW ' 
H^e /lave also the following estimates for the remainder terms 

\\R N \\ H » + \\S N \\ H s < MS N + 1 (N + l) 2 A N - 1 C^ +2 5 N + 1 e^ N + 1 '>^ t . 

Proof We have proven the H s estimates for all the terms of the expansion 
Uj, Vj, Tj, Qj. It is this easy to deduce, using (|73p . the estimate for the remainder 
terms. This comes from the inequality (1 < j < N — 1) 

\\D a (u j diu N - j )\\ < Cdl^HoollUjV-jllijIal+i + IK-HhII^Jv-jIIdo) 

(and subsequent inequalities), the Solobev embedding ||/||oo < II/II2 and the 
H s estimate for s = 2,3 for all the terms of the expansion, using also that the 
norm H s of the terms of the expansion in S 3 of order less than N is bounded 
by ClA^e^^K We thus deduce that 

A A \ (C AX\N-l„(N-l)y(k)t 

|| ^TjUh. < £ CA^C^e^ = CCoSe^ 1 -^ e . 

3=1 i=i 

When t < Tf = ^ In we obtain 1 - CoASe^W > 1-9, hence we deduce 

the estimate 

JV 



\\T N -l\\n.=\\^T j \\ H .<^6e«'-' 
3=1 

Moreover, one has 



N 

ir^-iiu^^ii^ii^-E^n^ 

3=2 



hence using 



N N 



Y.PWTjWL* < E^M*. ^ ^S 2 e 2 ^ 



3=2 3=2 



one obtains 



\\T N - 1|| L2 > SWT.mWe^ 3>^ e *r(*0t. 

One may thus consider C3 = We thus deduce that, for t < In ^^^j 1 ^ . 
we obtain 

||T JV -i|| La >^||r 1 (o)|U ie ^ t . 

Similar estimates hold for ||u H^a. 

Note that this proves that the first term of the expansion is the leading term of 
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the approximate total solution. 
For all what follows, we introduce 



I N+1 (t) = N ^A N - 1 C^+ 1 e^ N+1 ^ t . (77) 



5 Estimates of the (nonlinear) solution 

We constructed in the previous section a solution T N , u , Q N such that 
d t T N + u N WT N - k {x)u N T N = S N 

d t u N + u N .Vu N +T N p^W(p Q N ) = g + R N (78) 
divw^ = 

with the following properties for the remainder terms: 

| \pl d* n Rf\\ + | \pl 82« S N \\ < C n 5 N+1 I N+1 (t) (79) 

\\d^Rf || + \\d^S N \\ < C n 8 N+1 I N+1 (t) (80) 

i 

the constant C n depending on the Sobolev norm with weight p§ of the initial 
value of the normal mode solution and of the characteristic constants of the 
problem. 

_ i 

We deduced from this equality and the additional assumption kgp 2 bounded 
that we have identical estimates on R N and S : 

\\d^Rf\\ + \\d2«S N \\ < C n S N+1 I N+1 (t) (81) 

We study in this Section the global solution of the Euler system (|52l to obtain 
Sobolev estimates on the difference between the approximate solution and the 
full solution. Let T d = T - T N , u d = u ~ u N ,Q d = Q — Q N . We have the 
following system of equations: 



r d t T d + u N VT d + u d VT N = k a (uT d + u d T N ) - S N 

I Po(d t u d + u d Vu + u N Vu d ) + TV(poQ) - T N V( Po Q N ) = - Po R N (82) 
[ divu'' = 0. 

Before stating the results on the difference quantities according to the sys- 
tem, we use the properties of T N — 1, u N , Q N : 

Lemma 10 Let t e [0, T|]. 

For all a, there exists a constant C(|a|) such that 

\\D«(u d .Vu N )\\ < C(|a|)||tf d || w ^eT«* 

\\D a (T*(yQ N + koQ N ))\\ < C(\a\)\\T d \\ H ^e^ 

\\D°(k u d (T N - 1))|| < C(\a\)\\H% a{ ^fe^ 

\\D°(T N - l)(VQ d + feoQ^DH < C(|a|)||Q d || w+1 ^e^ fe )* 
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The proof of this Lemma comes from the fact that 

D a {fg N )=Y J C^fD^g N 

and we use the estimate H-D""^^^ < CHg^l^+iaH/ 3 !' as weu as ^ s 

result on any term of the form g N = ^2jLi$9j> wnere gj — uj,Vj,Tj,Qj to 

conclude for any term studied in the Lemma. Moreover, we use the Moser 
estimates to obtain 

\\D a (u.Vf) - u.VD a f\\ < C7(|]Vi?||oo]|V/|||«|_ x + IIV/HooPIh) 

hence, using u = u N + u d , one deduces 

and, similarity 

I ID" (T^Q^-T VD^Q**! | < C(5J(t) | IQ 11 ! 1 1„| +| | Vr* 1 ] loo 1 1 VQ d 1 1 1„| | VQ**! loo I 1 1 1 a , 

according to the equality D a {TVQ d ) - TVD a Q d = D a ((T - l)VQ d ) — (T — 
l)VD a Q d . 

We shall also use the following estimates 

||^(^.V^)||<C(|a|)||u d || 4 ||^||| Q | +1 (83) 

\\D a (u d .Vu d ) - tZ'.vWll < C(H)||^|| 4 ||^||| a |. (84) 

These equalities come respectively from (|74p and (|75|) . 

Introduce in what follows V = u d .Vu N + u N .Vu d , W = V+T d p~ 1 V(p a Q N ). 
We have the estimates 

\\D a V\\<M H I(t)5\\u d \\ H+ i 

\\D a W\\ < M lal I(t)d(\\v?\\ [al+1 + ||T d ||| Q |),Va 

\\D a (T d VQ d ) - T d VD a Q d \\ < C(||T d || 3 ||Q d || 2 + \\T%\\Q%) for |a| - 2 
\\T d p-h( P0 Q N )\\<5I(t)\\T d \\ 

5.1 Estimates on the density 

The equation on the density yield 

d t T d + u.X7T d - k Q uT d = k u d T N - u d NT N - S N . 

Apply the operator D a and denote by = D a {u.\JT d ) - u.WD a T d . This 
equation rewrites 

d t D a T d + u.VD a T d + W\ - D a (k Q uT d ) + D a {u d .VT N ) - D a (k u d T N ) = 0. 
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We can decompose — D a (kouT d ) into two parts, the one with u N , the other 
one with u d , denoted respectively by W a and W„ ■ It is clear that 

\\W» - D a (k Q uT d ) + D a (u d .VT N ) - D a k u d T N )\\ < C8(\\ti d \\\ a \ + \\T d \\ lal ). 

it is also clear that, using Moser estimates, ||W a || < CdlVu^Hooll^Hiai + 
1 1^*1 1 |q| Halloo)- O ne is thus left with the inequality 



\D a T d \\ < \\W a \\ + \\D a S N \\+CSI(t)(\\u d \\ H + \\T d \\ lal ). 



— I 

We have thus the estimate 

i\\D«T d \\ <C(||Vu d || co ||T d || H + ||^||| ct |||T d || co ) + ^+ 1 / JV+1 M 
+CSI(t)(\\u d \\ lal + \\T d \\ H ). 

5.2 Estimates on the pressure 

We obtained the relations 

||V0 d || < Midl^.V^H + SI(t)(\\u% + \\T d \\)+S N+1 I N+1 (t)). 

£|a=i HV^Q d || < M 2 (E, a | =1 \\D°(u d .Vu d )\\ + 8I(t)[\\u% + \\T%] 
+S N + 1 I N+1 (t) + (1 + SI(t) + \\T%)(\\u d .Vu d \\ 
+SI(t)(\\T d \\ + \\u d \\ 1 )+S N + 1 I N+1 (t)) 

Using the fact that t <T S , one obtains 

||VQ d || < M 1 (\\u d .Vu d \\ + H^lli + ||T d || + 5 N+1 I N+1 (t)). 

E|a=i HV^ a Q d || < M 2 (£ w=1 \\D«(u d .Vu d )\\ + \\u% + \\T% 

+ (1 + \\T%)(\\u d .Vu d \\ + \\T d \\ + \\u d \\ 1 ) + S N + 1 I N+1 (t)) 

In what follows, we introduce 

G% = D a (T N VQ d )-T N VD a Q d +D a (T N Q d k e 1 )+D a (T d VQ N +k T d Q N e 1 )+D a R N , 

G a = D a (T d VQ d ) - T d VD a Q d + D a (T d Q d k Q e 1 ). 
The equation on D a u d is 

d t D a u d + D a (u d .Vu d ) + D a R N + G a +G% + TVD a Q d = 0. 

When one multiplies by VD a Q d , one uses the divergence free condition on D a u d 
to get the estimate 

h\WD a Q d \\ < \\D a (u d .Wu d )\\ + \\D a R N \\ + \\G a \\ + \\G%\\. 

We use 

\\G%\\ < (7(1 + t)^ +3 (\\Q d \\ H + ||T d || H + \\u d \\ lal+1 )5I(t) 
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and 

||G Q || < C(\\VT d \U\Q d \\ H + ||7i| H ||VQ d |U + ||T d || 00 ||Q« , || w ). 
Hence we obtain (and it is pertinent for \a\ > 2) 

||V^g d || < C'(\\D<*(u d .Vu d )\\ + \\D°R N \\) + C5I(t)(\\Q d \\\«\ + \\T d \\\ a \ + ||w d || H+ i) 
+C(\\T%\\Q d \\ H + \\T d \\ lal \\Q%). 

For \a\ — 2, we will obtain ||Q d ||3, which is important. 
We use the equality, for \a\ =2 

||G Q || = \\D a T d VQ d + D f3 T d X7D a - f3 Q d CP\\ 

0</3<a 

which leads to the inequality 

\\G a \\<D (\\T d \\4\Q% + \\T%\\Q%). 
Replacing this estimate in the inequality for a such that \a\ = 2, one gets 
\\D a VQ d \\ < d (I C^T rf -^^T rf )| H-l Is^CH^I U| I i^^^^^^yv^i C*)- 
Using the inequalities on ||Q d ||i and ||Q d ||2, one gets 

||Q d ||i < Midl^.V^H + H^llr + ||T d || + 5 N+1 I N+1 (t)) 

||Q d || 3 < M 3 (||^.V^|| 2 + \\u% + \\T% + (1 + ||T d ||I + \\T%)\\u d .Vu d \\ + (1 + ||T d || 3 )||^|| 
+ ||T d || 4 ||^|| + S N ^I N+1 (t)(l + \\T d \U + (1 + ||T d || 3 ) 2 )) 

We use then the inequalities 

\\VD a Q d \\ < C{\\D a {u d .Vu d )\\ + HD^H + ||Q d || H 

+l|r d || H (i + ||Q d || 3 ) + II^IIm+i + (1 + l|r d ||3)||g d || H ) 

from which one obtains 

||Q d |lM+i < M N+1 (||« d .W d || H + II^Hh+i + S N+1 I N+1 (t) 

+r d || H (i + ||Q%) + ||Q d || H (i + ||T%)). 

Note that we have the estimate 

H^.WHh <C\\u%\\u d \\ lal+1 . (85) 

hence 

||Q d || H+ i <M w+1 ((l + ||tZ rf ||3)||t? , |||a|+i+« JV+1 /Ar + i(t) 

+l|r d ||| Q |(i + ||Q d ||3) + (i + ||t«|| 3 )||q<*|| w 1 j 
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It is then enough to use a recurrence argument to control the norm of Q d in 
H s+1 using the control of the norm of Q d in H". 

For the control on u d , let us rewrite the equation on D a u d . We introduce 

V a = D a (u.Wu d ) - u.WD a u d , W a = D a {T.VQ d ) - T.WD a Q d . 
We have the estimates 

\\V a \\ < C(l + H^IWH^IIh, \\W a \\ < C{1 + ||T d || 3 )||Q d ||| Q |. 
Using the relation 

J TVD a Q d D a u d dxdy = - J D a Q d (V(T N - 1) + VT d )D a u d dxdy 

thanks to the divergence free condition, as well as 

u.VD a v d .D a u d dxd// = i) 

one obtains the estimate 

j t \\D a u d \\ < ll^ll + ll^ll + ll^^ll + ll^^oTQ^II + II^Q^IKl + IIT^IIs), 
hence the inequality 

J t \\D a u d \\ < C[(l + ||T d || 3 )||g d |U + (1 + ll« d ||3)||^||| Q | +S N+1 I N+1 (t)} (87) 

For | a | > 3, this inequality is an a priori inequality. We have to state the 
identical inequalities for |a| = 0, 1, 2. 
We have the following inequalities: 

j t \\u d \\ < C ((l + \\T%)\\Q d \\ +S N+1 I N+1 (t)) (88) 

because V a — W a = 0, 

^ 1 | < ^-(Ci-^l l^^l | 3 )l li^Ci^l Is) 1 1^| U^Ci^l U)l la^^^^ 1 ^^^! C*)) 

(89) 

and 

^ 1 1-^^ I | s < ^-CCi-f-l | 3 ) I I I^^Ci^l I 1 1 i^Cl^l la) I l^^l l^^^^^^^r^^i C*)) 

(90) 

We thus deduce an estimate of the form 

j t (\\T d \\l+\\u d \\l) < <C7Ci^| U) 4 ^| |^| | 3 )CI 12-^1 1^1 1^)^^^^ 1 ^^-^! C*)CI l^^l 11^1 

from which one deduces an estimate of the form 
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End of the proof We thus know that, for t < T 5 , we have d N+1 I N+1 (t) < 1 
hence an inequality of the form 

jH{t)<C{(l + (H(t)f)H(t) + l) 

where H(t) = {\\T d \\\ + \\u d \\l)h . 

As we have H(0) = 0, one deduces that 



H{t) ds 
o (l + ^ + l" 

The function H — > J Q H nrp^i+i * s a bijection from [0, +oo [ onto [0, Jq + °° ^ 1+ ^ s+1 [■ 
For > 1, one deduces Ct > f Q (1+f f/ )s+1 , hence for t < ± J Q = 
T\, one obtains H(t) < 1. The set of points t such that t > and iJ(i) < 1 is 
not empty. 

Once this set is not empty (and once we proved that the solution exists for a 
time Ti), we obtain 

Lemma 11 Let h be a function such that 

^<C(l + h(t)) 4 h(t) + C5 N+1 eS N+1 ^\ h(0) = 0. 

For S < 1 and (N + l)j(k) > 17C, denoting by T 5 = ^ In \, one has 

Vte[0,T*],h(t)<5 N+1 e( N+1 ^ t . 

Proof The inequality we start with is 

j t h{t) < C{1 + h(t)fh(t) + C8 N+1 e {N+1 ^ t . 

We consider N such that (N + l)-f(k) > 17C. We study the interval where 
h(t) € [0,1], knowing that h(0) = 0. Consider to the first time (if it exists) 
where h(t ) = 1. If it does not exist, then h(t) < 1 for t e [0, Tq] and we have, 
for all t G [0, T 5 ] the inequality 

h'(t) < 16Ch(t)+CS N+1 e^ N+1 ^ t . 
from which one deduces 

h(f) < _ (W+l) 7 (fe)t WV+l (JV+l) 7 (fe)t 

m S (JV + 1)A-16C 
hence h(T$) < 1. 

If to exists, we have, for all t € [0,io]> the inequality 

±(h(t)e- wct ) < -C(l - h(t))h(t)R(h(t))e- 16Ct +CS N + 1 e^ N+1 ^ t - 16Ct 
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where R(x) = (1 + x) 6 + 2(1 + x) + 4(1 + x) + 8, from which one deduces that 

[ 0) S (N + l) 7 (fc) - 16C < 

hence /i(io) < L contradiction. 

We thus deduce that h(t) < 1 for t G [0,T 5 ], hence 

^)<^ +1 e™^fe[o,r j ]. 

Lemma [TT] is proven. 

We have thus the inequalities 

„2 7 (fe)t 

INI > IKH - I Nil > s\\mo)\\ - c o AS \_ Co6Ae7(k)t ~ ^ +le(N+lh(fc) - 



Choose i = T? = ^ In We have 

I Nil >fc^* [|K(o)| | -Co 



1 - 



We thus check that there exists €q < | such that 8 < Eq implies [| 
Coj^g ~d N ] > ||Ni(0)||. Hence for i < ^ln^, one has 

IN(t)||>^INi(o)||^( fe )*. (9i) 

In particular 

We proved Theorem O 

It is then clear that, for < t < Tf, this term is smaller than #,as small as one 

wants, hence the inequality on T d , u d . 

As T = T N + T d , u = u N + u d , one obtains 

llT-llU^llT^-rt 

which imply the result. 
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